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Abstract 

This work is devoted to giving a geometric framework for describing higher-order non- 
autonomous mechanical systems. The starting point is to extend the Lagrangian-Hamihonian 
unified formalism of Skinner and Rusk for these kinds of systems, generalizing previous devel- 
opments for higher-order autonomous mechanical systems and first-order non-autonomous 
mechanical systems. Then, we use this unified formulation to derive the standard La- 
grangian and Hamiltonian formalisms, including the Legendre-Ostrogradsky map and the 
Euler-Lagrange and the Hamilton equations, both for regular and singular systems. As ap- 
plications of our model, two examples of regular and singular physical systems are studied. 
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1 Introduction 

Higher-order dynamical systems play a relevant role in certain branches of theoretical physics, 
applied mathematics and numerical analysis. In particular, they appear in theoretical physics, 
in the mathematical description of relativistic particles with spin, string theories, Hilbert's 
Lagrangian for gravitation, Podolsky's generalization of electromagnetism and others El [71 El 
[271 [321 [Ml [371 [38l Uni HTj , as well as in some problems of fluid mechanics and classical physics 
(see, for instance, the example in Section [6. II taken from [21 [2S]), and in numerical models arising 
from the discretization of first-order dynamical systems that preserve their inherent geometric 
structures [21]. In these kinds of systems, the dynamics have explicit dependence on accelerations 
or higher-order derivatives of the generalized coordinates of position. 

In recent years, much works has been devoted to the development of geometric formalisms for 
higher-order mechanics and field theory (see, for instance, [Il[21[IIl|12l[l6l|22l[26l[28ll2l[3ll|33l 
144]). These formulations use higher-order tangent and jet bundles as the main tool. In particular, 
in a recent paper [39j a new geometric formulation has been proposed, which is an extension to 
higher-order autonomous mechanical systems of the formalism proposed by R. Skinner and R. 
Rusk in his seminal paper [l5]. This formulation compresses the Lagrangian and Hamiltonian 
formalisms into a single one, originally developed for first-order autonomous mechanical systems 
and later generalized to non-autonomous systems (5] [T3], control systems [1], and first-order 
classical field theories (see [41] and references therein). Nevertheless, to our knowledge, there 
is neither a complete geometrical description of the Lagrangian and Hamiltonian formalisms 
(partial studies on this subject can be found in [151 [SOI [231 UHl EO])) nor of the Skinner- Rusk 
unified formalism for non-autonomous higher-order mechanical systems. 

The aim of this work is to fill this gap. In order to do this, we first develop the Lagrangian- 
Hamiltonian unified formalism of Skinner- Rusk for higher-order non- autonomous mechanical sys- 
tems, studying in particular how this formulation enables us to obtain the generalized Legendre- 
Ostrogradsky map connecting the Lagrangian and Hamiltonian formalisms, as well as the Euler- 
Lagrange and the Hamilton equations of motion. Thus, starting from this unified framework, we 
obtain a geometric description for the Lagrangian and Hamiltonian formalisms for higher-order 
non-autonomous mechanical systems. This study is conducted both for regular and singular 
dynamical systems. Our analysis is performed by using higher-order jet bundles, since we wish 
this work to serve as a model to develop an unambiguous framework for higher-order classical 
field theories that complete previous approaches in this way [101 146] . 

The paper is organized as follows: in Section [2l we review the geometric structures needed 
to develop the formalism, such as the higher-order jet bundles, the total derivatives and higher- 
order semisprays. Section [3] is devoted to the geometric formulation of the Skinner-Rusk unified 
formalism for higher-order non-autonomous mechanical systems, including the description of 
the dynamical equations using sections and vector fields. In Sections [H and [5] we recover the 
standard Lagrangian and Hamiltonian formalisms, presenting a complete description of both for 
regular and singular systems. Finally, in Section [6l two examples are analyzed; the first is a 
regular system which models the shape of a deformed elastic cylindrical beam with fixed ends 
and has applications in Statics and other branches of classical physics [H [25] ; the second is a 
modification of a singular system describing a relativistic particle [36l [351 El IMl [39] , which in our 
case is subjected to a generic time-dependent potential. The paper concludes in Section [7] with 
a summary of results and future research, and an appendix in Section [A] where the particular 
situation of higher-order trivial bundles is briefly analyzed. 

All the manifolds are real, second countable and C°°. The maps and the structures are 
assumed to be C°°. Sum over repeated indices is understood. 
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2 Geometric structures of higher-order jet bundles over M 

2.1 Higher-order jet bundles over M 

(See [221113] for details). 

Let E M be a bundle (dimi? = n+l), and let rj € be the canonical volume form in 

R. If k € N, the kth order jet bundle of the projection vr, J^tt, is the {{k + l)n + l)-dimensional 
manifold of the A: -jets of sections (j) G r(7r). A point in j'^vr is denoted by where (f) S r(7r) is 
any representative of the equivalence class. We have the following natural projections: if r ^ A;, 

vr^: J'^TT , TT^: J^tt — > E 

fcj) I — > fcj) , fcl) I — > 4) , 

Notice that ttq = vr'^, where J^tt is canonically identified with E, and 7r| = Idjk^. Furthermore, 
we denote tt^ = -k o -k^ : J^ir — )- M. 

Local coordinates in J^tt are constructed as follows: let t be the global coordinate in R such 
that r] = dt, and (t, g"^), (1 ^ A ^ n), local coordinates in E adapted to the bundle structure. 
Let (f) G r(7r) such that (f) = {t, (f>^). Then, local coordinates in j'^vr are (t, q"^, qf^, . . . , q^), with 

Usually we write q^ instead of q^, and so the local coordinates in J^ir are written {t, q^, qf, . . . , q^) 
Using these coordinates, the local expression of the natural projections are 

7r^(i, Qo. 9^ • • • , Qk) = it, qo, 9^ • • • , Qr) , Qo^ 9^ • • • > Qk) = (t^ Qo) ■ 

If i;^ G r(7r) is a section of vr, we denote by the canonical lifting of (f> to j'^vr, that is, the 
map j'^fj)' ^ ~^ J^TT, which is a section of the projection vf^. 

Remark: We use the same notation for points of j'^vr and liftings of sections to j'^vr, since 
giving a point in J^vr is equivalent to giving the lifting to J*^7r of a section of vr (see for 
details) . 

2.2 Total time derivative 
(See |l3] for details). 

Definition 1 Let E M 6e a bundle, to G M, (/> G r(vr) and u G T^^R. The kth. holonomic 
lift of u by 4) is defined as ((jV)*(ii), jf/V) G (7r^+^)*Tj'=vr , where jl+^(t> = {f^^(p){to). 



d 

In local coordinates, if u is given by ti = Uqtt 

at 



, the kth holonomic lift of u is given by 

to 



k 



The vector space (vr^"''^)*(Tj'^vr) .fe+i has a canonical splitting as a direct sum, as follows: 
(vr^i)*(Tj'=^) = (vr,^+i)*(F(vf'=)) e (jV)*Ti„R, 



(1) 
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where (j i?i))*T(^]R denotes the set of kth. holonomic Hfts of tangent vectors in T^^M by 0. As 



a consequence, the vector bundle (^^+i)*Tj'^7r " '"'^ vr may be written as the direct 
sum of two subbundles: 

where H{-k^'^^) is the union of the fibres (j^(/))*T(M, for t G M. 

Now, if X(7r^'''^) denotes the module of vector fields along the projection vr?"*"^, the submodule 



corresponding to sections of i'^'l~^^)*Tjk^ _ is denoted by X^{it^'^^), and the submodule 



corresponding to sections of (vr^^^)*rjfc^ ^ ^^^^ is denoted by ^''(vr^'''^). The splitting for the 
bundles given above induces the following canonical splitting for the module ^(vr^'''^): 

An element of the submodule X^{Tr^'^^) is called a total derivative. 

Definition 2 Given a vector field X € jC(M), a section (p G r(7r) and a point to G M, the ki\i 
holonomic lift of X by cp, X^ = j'^X £ X^{tt'^'^^), is defined as 

Hence, every vector field X G jt(R) defines a total derivative given by its holonomic lift. 

Alternatively, we have the following characterization of X^ as a derivation: for every / S 
C°°(j'=7r) we have 

where dxk is the derivation associated to X^ and dx is the derivation corresponding to X. 

d 

In local coordinates, if X G X(M) is given by X = Xq-—, then, bearing in mind the local 

ot 

expression of the kth holonomic lift for tangent vectors ([T]), the kth holonomic lift of X is 

Finally, the total time derivative is the fcth holonomic lift of the coordinate vector field 
d/dt £ X(M), which is denoted by dr G X{7r^~^^), and whose local expression is 

d A A d 



A 



1=0 * 

Remark: The usual notation for the total time derivative is d/dt, as seen in |43j . while the 
notation dx is usually reserved for the same operator in the autonomous case. Nevertheless, in 
this paper we use the same notation for both operators, and the one that is considered will be 
understood from the context. 
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2.3 Higher-order semisprays. Holonomic sections 

Now we generalize the concept of semispray introduced in [23 to the time-dependent case. 

Definition 3 A section if) G r(7f'^) is holonomic of type r, 1 ^ r ^ A;, if j^^^'^^cj) = vr|_,^_(_^ o i/;, 
where (f> = tt^ o ip ^ r(7r); that is, the section tj) is the lifting of a section of t: up to J'^~'^^^'it. 




In particular, a section ip is holonomic of type 1 if, with (p = tt^ o ij; , then j^cp = ip; that is, 
tp is the canonical k-jet lifting of a section (p S r(7r). Throughout this paper, sections that are 
holonomic of type 1 are simply called holonomic. 

Definition 4 A vector field X G X( J'^vr) is a semispray of type r, 1 ^ r ^ k, if every integral 
section tp of X is holonomic of type r. 



The local expression of a holonomic section of type r, tp ^ T{J^it), is 

Thus, the local expression of a semispray of type r is 

d A d d 



X 



+ Q2 



+ 



+ Qk-r+l 



+ ^k-r+1 



d 



Prom the local expression, it is clear that every holonomic section of type r is also holonomic 
of type s, for s ^ r. The same remark is true for semisprays. 

We observe that, from the definition, semisprays of type 1 in J^ir are the analogue to the 
holonomic vector fields in first-order mechanics; that is, they are the vector fields whose integral 
sections (curves) are the canonical liftings to J^ir of sections (curves) on the basis. Their local 
expressions are 



+ 92 



d 



+ • • • + C 



d 



Xi 



d 



If X G X( J'^vr) is a semispray of type r, a section (p E r(vr) is said to be a path or solution of 
X if j''(p is an integral curve of X; that is, j'^cp = X oj^cp, where j'^cp denotes the canonical lifting 
of j^(p from J^TT to T(j'^7r). Then, in coordinates, (p verifies the following system of differential 
equations of order k + 1: 



^k-r+1 



dt ' ■ ■ ■ ' dt* 



Xk ( 
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3 Skinner-Rusk unified formalism 



3.1 Unified phase space. Geometric and dynamical structures 



Consider the configuration bundle vr : — )■ M, where E is an (n-l-l)-dimensional smooth manifold. 
Let C G fi^{J^TT) be a kth order Lagrangian density, that is, a vf'^-semibasic 1-form. Thus, we 
can write C as C = L - [Tt^)*r} = Ldt , where L G C°°( j'^vr) is the kth-order Lagrangian function. 

According to [SIIMIE^, we consider the following bundles: 



J vr Xjfc-ij 



J vr 



(the fiber product of the above bundles), where J^~^tj:* = T*( j'''^^7r)/(7f'^^^)*T*M. The bundles 
W and Wr are called the higher-order extended jet-momentum bundle and the higher-order 
restricted jet-momentum bundle, respectively. 

Comment: The reason for taking these bundles is in order to recover the Lagrangian and 
Hamiltonian formalisms from this unified framework, and as we see in Sections H] and \5\ those 
formalisms take place in the bundles J'^'^~^it and J^~^tt*. 



These bundles are endowed with the canonical projections 



J vr 



— >■ J TT 



PR - W 
pi- Wr - 



In addition, the natural quotient map p: T*(j'^ ^vr) — )• J*^ ^vr* induces a natural projection 
(that is, a surjective submersion) py\; : W — ?■ Wr- Thus, we have the following diagram 




T*(J 













where ttj^-i^: T*(J vr) — )■ J vr is the canonical submersion and vr^fc.i^: J' 
is the map satisfying njk-i^ = tTj^^i^ o p. 



If {U; t, Qq) is a local chart of coordinates in E, we denote by ((vr 
and ((tt^.-i, ovr'=-i)*(C/);t,g^, 



2k-l\ 



'{uy,t,q^ 



J 



' ^2k~l) 

iQk-i^P^PA' ■ ■ ■ ^Pa'^) the induced local charts in J^^~^7r 
and T*(j'^~^7r), respectively. Thus {t,qQ,...,q^_-^^,p^,...,p^^) are the natural coordinates 



in J vr*, and the coordinates in W and Wr are (t, Qq, . . . , q^_i. 



'<l2k-vP^P% 



fc-l^ 

'Pa > 
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and {t, Qq, . . . , q^_i, q^, • • • , Q2k-i^PA^ ■ ■ ■ respectively. Note that dim W = S/cn + 2 and 

dim(>Vr) = S/cn + 1. 

The bundle W is endowed with some canonical geometric structures. The first one is: 

Definition 5 Let € i7^(T*( J^^^vr)) he the tautological 1-form, and Qk-i = ~d&k-i S 

i7^(T*(J*'~^7r)) the canonical symplectic 2-form on T*[J^~^t:). We define the higher-order 
unified canonical forms as 

e = p*efc_i G ; = Ps^fc-i e . (3) 

Bearing in mind that the local expressions for the canonical forms on T*(j'^^^7r) are 

Qk~i=p\dqt +pdt ■ VLk-i=dqt ^dp\-dp^dt, (4) 
the above forms can be written locally as 
= P*2{p\Mt+P<^t) = p\dqt+pdt ■ n = p*{dqfAdp\-dpAdt) = dqfAdp\-dpAdt (5) 

Notice that from the local expressions ([5]) we have 



Thus, is a presymplectic form in W. 

The second canonical structure in W is the following: 

Definition 6 The higher-order coupling 1-form in W is the p^-semibasic 1-form C G i7^(W) 
defined as follows: for every w = {y,aq) G W (that is, aq G T*{J^~^tt), where q = '/r^^"^^(y) is 
the projection of y to J^'^tt) and u G T^W, then 

{C{w) I u) = {aq I (T^(/-Vom))W>, (7) 
where 4> G r(7r) is any representative of y (that is, p^~^(f) = y). 

C being a pK-semibasic form, there exists C G C°°(yV) such that C = Cp^r] = Cdt. An easy 
computation in coordinates gives the following local expression for the coupling 1-form: 

C = (p + pkii)dt. (8) 

We denote C = (vrf'"^ o p^)*C g f2^{W). As the La grangian density is a vr -semibasic 
form, we have that £ is a /3iR-semibasic 1-form, and thus we can write C = Lp^ij = Ldt, where 
L = (vr^'^"^ o pi)*L G C°°(W) is the pull-back of the Lagrangian function associated with C. 
Then, we define a Hamiltonian submanifold 



Wo 



ju- G W: C{w) =C(tf)} Aw. 



C and £ being /OK-semibasic 1-forms, the submanifold Wo is defined by the constraint C — L = 0. 
In natural coordinates, bearing in mind the local expression ([8]) of C, the constraint function is 

C-L=p + p\qf^^-L = 0. 
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We have the following natural projections in Wo- 

Local coordinates in Wq are {t, Qq, . . . , q^-i^Qk^ ■ ■ ■ ■> 92fe-i'^'A' • • • jPa"^)' ^^"^ local expres- 
sions of the above maps are 

Pl{t, Qf, (if, Pa) = (t, <lf, Qf) ; P2(*> <lt^4M = (t, Qf, L - p^qf+iM 
P2{t, Qf, if, Pa) = {t, Qf,PA) ; jo{t, qf, qf,p\) = qf, qf, L - p\qf+i,p\) ■ 

Proposition 1 The submanifold Wo ^ W is 1-codimensional, iJ,y\> -transverse and diffeomor- 
phic to ■ 

(Proof) Wo is obviously 1-codimensional, since it is defined by one constraint function. 

To see that Wo is diffeomorphic to Wr , we show that the smooth map /Uyy o : Wq — )■ Wr is 
one-to-one. First, for every (y, a) G Wo, we have L{'K^~^{y)) = L{y,a) = C{y,a), and 

{m°3o){y,Oi) = IJ'W{y,a) = {y,n{a)) = (y, [a]) . 

First, /i>v ° jo is injective; in fact, let (yi,ai), (^2,0:2) € Wo, then we wish to prove that 

(mw o Jo)(yi,ai) = (m oio)(y2,a2) <^ (yi,ai) = (^2,0:2) <^yi=y2 and ai = 02 . 

Now, using the previous expression for (//yy ° jo){y,Oi), we have 

[fiy^o jo){yi,ai) = (mw o Jo)(y2, 0:2) <^ (yi, [ai]) = {m, [q!2]) <^yi = y2 and [ai] = [0:2] . 

Hence, by definition of Wo, we have L(7r^''"^(yi)) = L(7r^^"^(y2)) = C{yi,ai) = C{y2,a2). 
Locally, from the third equality we obtain 

p{ai) + p\{ai)qf_^i{yi) =^(02) + PA{a2)qf+i{y2) , 

but [ai] = [02] =^ =^^(["1]) =Pa([«2]) =Pa("2)- Thenp(ai) =p{a2), and ai = 02. 

Furthermore, /iyv ojo is surjective. In fact, given (y, [a]) G Wr, we wish to find (y,/?) G Jo(Wo) 
such that [j3] = [a]. It suffices to take [P] such that, in local coordinates of W, 

pm=p'Am) , piP) = L{4'-\y)) - PAi[^])<itiiy) ■ 

This P exists as a consequence of the definition of Wq. Now, since /xvv°Jo is a one-to-one submer- 
sion, then, by equality on the dimensions of Wo and Wr, it is a one-to-one local diffeomorphism, 
and thus a global diffeomorphism. 

Finally, in order to prove that Wq is /xvv;-transversal, it is necessary to check if L(Y){^) = 
y(^) 7^ 0, for every Y G ker//w^ and every constraint function ^ defining Wo- Since Wo is 
defined by the constraint function C — L = and ker //>v^, = {d/dp}, we have 

^ — - ^ ' J ^A 



then Wo is //yy-transversal. 
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As a consequence of this last result, in the following we consider the diagram: 

W. 




T*(J 













As a consequence of Proposition [U the submanifold Wo induces a section /i G r(^H'), that 
is, a map /i: Wr — W. This section is specified by giving the local Hamiltonian function 



H = -L + p\q:, 



i+l ) 



(9) 



that is, h{t , qf , , p\) = {t,qf,q^,—H,pl^). The section h is called a Hamiltonian section of 
fj^yv, or a Hamiltonian fj.\A^;-section. 

Next, we can define the forms 

with local expressions 

So = p\dqf + {L - p\qf_,,)dt ; = dqf A dp\ + d{p\qf^, - L) A dt , (10) 

and we have the presymplectic Hamiltonian systems (Wo,f^o) and (yVr,^r), with 0^ = 
Finally, it is necessary to introduce the following concepts: 

Definition 7 A section ipo G is holonomic of type r in Wo, l^r^2k — l, if the section 

Pi o ipg £ r(7f^'^~"^) is holonomic of type r in J'^^~^it. 

Definition 8 A vector field Xq € X(Wo) is said to he a semispray of type r in Wo if every 
integral section ipo of Xq is holonomic of type r in Wo- 



The local expression of a semispray of type r in Wo is 



Xo 



j=0 ^« i=2k-r ^« i=0 

and, in particular, for a semispray of type 1 in Wo we have 

2k-2 ^ o fc-1 



d 



j=0 



^'i+l'fTA" -^2k-l o A 



d 

^dp\ 
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3.2 Dynamical equations 

The dynamical equations for non-autonomous dynamical systems in general can be geometrically 
written in several equivalent ways, using sections (curves) which are the dynamical trajectories, 
or vector fields whose integral curves are the dynamical trajectories. In this section we explore 
both of these ways, and prove their equivalence. 



3.2.1 Dynamical equations for sections 

The Lagrangian-Hamiltonian problem for sections associated with the system {Wo,^o) consists 
in finding sections ipo G r(/9j|) (that is, curves V'o : K — )• Wo) characterized by the condition 

i>l i{Y)no = 0, for every Y G X(Wo) . (11) 

In natural coordinates, let Y G X(VVo) be a generic vector field given by 

dt ^ dqf dqf ^ dp\ 

bearing in mind the coordinate expression (jlOp of i^o? the contraction is 

+ ff [^P\ - §^dt + p^'dt^ + F,^ (^p'X~' - ^) dt + G\ {qt.dt - dqf) . 
Thus, taking the pull-back by the section V'o = {t,qf'{t),qf{t),p\{t)), we obtain 
C i{Y)no = f (^^qf - qt,p\ - P\it^ dt + ff - ^ dt 

+ ft{p\-^+ p'a"] dt + [pI'' - ^) + i^ti - it) dt . 



Finally, requiring this last expression to vanish and bearing in mind that the equation must 
hold for every vector field Y G X(Wo) (that is, it must hold for every function f,ff',Fj^,G\ G 
C°°{Wo)) we obtain the following system of equations 



dL 

qt-qf+iPA-PAQf+i = o (13) 



dqf 

p'a = ^ (14) 
dqf 

Br 

if = lAi ■ (17) 

It is easy to check that equation (|13p is redundant, since it is a consequence of the others. 
Equations , and p!7|) are differential equations whose solutions are the functions defining 
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the section ipo- In fact, equations ([HD, (fTC|) give the higher-order Euler-Lagrange equations, as 
we see at the end of this Section and in Section [H In addition, observe that equations (jl6p do 
not involve any derivative of ipo- they are pointwise algebraic conditions. These equations arise 
from the /52-vertical part of the vector fields Y. Moreover, we have the following result: 

Lemma llfYe X^(''2)(V\;q)^ then i(y)Oo is p^-semibasic. 

(Proof) A direct calculation in coordinates leads to this result. Bearing in mind that a local 
basis for the vertical vector fields is given by ([6]) and the local expression (jlOp of r^o, we have 



d 



no 



PA'-Q^jdt, for J = k; 



dL 



= 0-dt, ior j = k + l,...,2k-l. 

Thus, in both cases we obtain a /o^-semibasic form. 

As a consequence of this result, we can define the submanifold 

Wc = {w(^Wo: {i{Y)no){w) = for every Y G ^^('^^^(Wo)} A Wo , 



where every section ■0o solution of equation must take values. It is called the first constraint 
submanifold of the Hamiltonian presymplectic system (Wo,f^o)- 

Locally, Wc is defined in Wo by the constraints — dL/dq^ = 0, as we have seen in (fTHl) 
and in the proof of the previous Lemma. In combination with equations (jlSp . we have: 



Proposition 2 Wc contains a submanifold Wi ^ Wc which can be identified as the graph of a 
map TL: J^~^tx* defined locally by 

(Proof) As Wc is defined locally by the constraints — dL/dq^ = 0, it suffices to prove 
that these constraints give rise to the functions defining the map given above, and thus to the 
submanifold Wi. We do this in coordinates. 

Taking into account that dT(p\) = p\ along sections, the constraint function defining Wc, 
in combination with equations (jlSp give rise to the following constraint functions 



dL 



Pa-[^- Mpa) 1 = Pa 
L 



Pa- ^^- dT(p\) ]=p'a- > I .f^ 1 = , 
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Therefore, these constraints define a submanifold Wi ^ Wc and we may consider that this Wi 
is the graph of a map TC: j'^^^^ir — t- J^^^tt* given by 



dL 



r+i 



Bearing in mind that the submanifold Wo W is defined locally by the constraint function 
P + p\(lt+i — L = 0, and that Wi is a submanifold of Wc, and thus a sumbanifold of Wo, from 
the above Proposition we can state the following result, which is a straightforward consequence 
of the previous result: 

Corollary 1 Wi is the graph of a map TL: j'^^-^Ti T*(j'=~i7r) defined locally by 

T£,*t = t , TC*qf = q^ , 



Remark: The submanifold Wi can be obtained from Wc using a constraint algorithm. 
Hence, Wi acts as the initial phase space of the system. 

The maps J-'C and J-'C are called the extended Legendre-Ostrogradsky map and the restricted 
Legendre-Ostrogradsky map associated to the Lagrangian density C, respectively. A justification 
of this terminology is given in Section [5j Now we can give the following definition: 



Definition 9 A Lagrangian density C G f2^{J^7r) zs regular if the restricted Legendre-Ostrogradsky 
map T C is a local diffeomorphism. Lf the map J-C is a global diffeomorphism, then C is said to 
be hyper regular. 



Computing in natural coordinates the local expression of the tangent map to J^£, the regu- 
larity condition for C is equivalent to 



d'^L 

det ( q^bq^a ) (y) ^ ^^'^^y y ^ -^^^ ■ 



Equivalently, if we denote ff^^ = TL*p^^^ ^ then the Lagrangian density L is regular if, and 
only if, the set (t, g/^,j5^), O^i^/c — 1, isa set of local coordinates in J^^^^vr. The local 
functions are called the Jacohi-Ostrogradsky momentum coordinates, and they satisfy that 

.1 _ dL 
dqi. 

which are exactly the relations given by (jlSp . taking into account that dT = d/dt along sections. 



Pa' = 7^- MP a) , (18) 



Notice that equations ([HI), p^ . and pT|) do not allow us to determinate the functions g^, 
A; ^ j ^ 2fe — 1, of the section tpo. Thus, in the general case, we need an additional condition 
when stating the problem, which is the holonomy condition for the section V'o- Therefore, the 
Lagrangian-Hamiltonian problem must be reformulated as follows: 

The Lagrangian-Hamiltonian problem consists in finding holonomic sections ipo £ ^(/'r) charac- 
terized by the equation pT]) . 



Remarks: 
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In fact, the functions q^, k ^ j ^ 2k — 1, are determined by the equations (|14|) and (|15p . 

bearing in mind that the section ■0o must He in the submanifold Wi = graph(J^£). It is 
easy to see that, by replacing the local expression of the extended Legendre-Ostrogradsky 
map in the equations (jl4p and (|15p . these equations lead to the Euler-Lagrange equations 
and to the remaining {k — l)n equations that give the full holonomy condition: 



dL 



dq^dq^ 
d dL 



i=k 



' )(• 



li+l 



(A; ^ j ^ 2A; - 2) 



■>Po 



dt dq^ 



+ 



d? dL 



dt2 dq^ 



+ 



+ (-1) 



d'^ dL 



dt^ dq^ 



i'o 



where the terms in brackets (• • • ) contain terms involving partial derivatives of the La- 
grangian function and iterated total time derivatives, and the first sum (for j = k) is 
empty. However, observe that these equations may or may not be compatible, and a suffi- 
cient condition to ensure compatibility is the regularity of the Lagrangian density. Thus, 
for singular Lagrangian densities, the holonomy condition for the section ipo is required. 

The requirement of the section ipo to be holonomic is a relevant difference from the first- 
order case, where the holonomy condition is deduced straightforwardly from the dynamical 
equations when written in local coordinates. Nevertheless, in the higher-order case, the 
equations allow us to recover only the holonomy of type k, as seen in (I17p . and the highest- 
order holonomy condition can only be recovered from the equations if the Lagrangian 
density is regular. Hence, this condition is required "ad hoc". 

The regularity of the Lagrangian density has no relevant role at first sight. However, as we 
have seen in the first remark, equations (jl4p and (jlSp give the higher-order Euler-Lagrange 
equations, which have a unique solution if the Lagrangian density is regular. For singular 
Lagrangians, these equations may give rise to new constraints, and a constraint algorithm 
should be used for finding a submanifold where the equations can be solved. 



3.2.2 Dynamical equations for vector fields 



The Lagrangian-Hamiltonian problem for vector fields associated with the system (Wo,^o) 
consists in finding vector fields Xq G X{yVo) such that 



(19) 



According to [13j and [i8\, we have: 



Proposition 3 A solution Xo G X(Wo) to equation (fT9|) exists only on the points of the sub- 
manifold Sc defined by 



5c = ju; G Wo: {i{Z)dH){w) = for every Z G ker J^} A Wo , 



(20) 



We have the following result: 



Proposition 4 The submanifold Sc ^ Wo contains a submanifold Si ^ Sc which is the graph 
of the extended Legendre-Ostrogradsky map; that is, Si = graph and hence Si = Wi. 
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(Proof) As Sc is defined by (j20p . it suffices to prove that the constraints defining Sc give 
rise to the constraint functions defining the graph of the extended Legendre-Ostrogradsky map 
associated to £. We do this calculation in coordinates. Taking the local expression ([9]) of the 
local Hamiltonian function H G C°°(Wo), we have 

i=0 1=0 '^i 

and using the local basis of ker $7 given in ([6]), we obtain that the equations defining the sub- 
manifold Sc are 

dL 

i{Z)dH = p^-^ J = , for every 1 ^ ^ ^ n . 

Note that these expressions relate the momentum coordinates with the Jacobi-Ostrogradsky 
functions j5^~^ = dL/dq^, and so we obtain the last group of equations of the restricted 
Legendre-Ostrogradsky map. Now, using the same argument as in the proof of Proposition 
[2] and the relations p8|) for the momenta, we can consider that Sc contains a submanifold Si 
which is the graph of a map 

{t,qf,qf) ^ {t,qf,p\) 
which we identify with the restricted Legendre-Ostrogradsky map by making the identification 

r — 1 — 1 

Pa =Pa ■ 

Finally, taking into account that Si is also a submanifold of Wo, which is defined by the 
constraint p + p\q{^i — L = 0, we have as a direct consequence that Si is the graph of the 
extended Legendre-Ostrogradsky map TC, and hence Si = Wi. ■ 

We denote by Xw^CW) the set of vector fields in Wo at support on Wi. Hence, we look for 
vector fields Xq G Xwj(Wo) which are solutions to equations ([1^]) at support on Wi; that is 

i(Xo)17olwi = ; i(Xo)(p^)*77l>v, = 1 • (21) 

In natural coordinates, let Xq G X(Wo) be a generic vector field given locally by p^ . Thus, 
from (jl9p we obtain the following system of {2k + l)n -|- 2 equations 

cA ^-^ , fA i-1 dL \ ^ ui dL \ j A 



= fiti , (23) 

Pa' = fMp\) , (24) 




\dqf J 

(25) 
(26) 

where O^i^/c — lin ()23p , and l^i^A; — lin (j24p . By a simple calculation one can see that 
equation ()22p is redundant, since it is a combination of the others. Therefore 

d A d d dL d ,,4.9 , „, 
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Remark: In a more general situation, the second equation in (|19p is written i{Xo)p^r] ^ 0, 
that is, a p^-transversal condition for the vector field Xo. In local coordinates, this replaces 
equation ([2^ by / 7^ 0, thus giving the vector field 




where / G C°°(Wo) is any non- vanishing function. This gives a whole family of vector field 
solutions to the dynamical equations, and taking a particular constant value for / just fixes a 
specific vector field in this family. From a physical viewpoint, taking a particular value for / is 
just fixing the gauge. 



Observe that equations (j26p are just a compatibility condition for the vector field Xo, which, 
together with the relations (|18p for the momenta, state that vector field X^ solutions to equations 
([T9|) exist only at support on the submanifold defined by the graph of the extended Legendre- 
Ostrogradsky map. Thus, we recover, in coordinates, the result stated in Propositions [3] and [H 
Furthermore, equations (|23p show that Xq is a semispray of type k in Wq. 




The component functions F^, k ^ j ^ 2k — 1, are undetermined. Nevertheless, recall that 
Xo is a vector field that must be tangent to the submanifold Wi. Thus, it is necessary to impose 
that L(-^o)CIwi ~ every constraint function ^ defining Wi. Locally, this is equivalent to 
imposing ^o(C)lwi ~ Hence, taking into account Prop. HI these conditions lead to 

I d A d ^A d dL d , , 



( d A 9 ^A d dL d , , ,• 

(d A d ^A d dL d , , ,• 

(observe that we do not need to check L(^o)(p — FC p) = 0, since this is the constraint defining 
the submanifold Wo ^ W, and Xo is a vector field already defined in Wo) and, from here, we 
obtain the following kn equations 

-,B R \ d^L R \ , / d^L 




{Fk+i - Qk+2) jr^mni - i^k - 1k+i) dr 



5« ydQ^dql 
d'L ^^„R R 



dq^dq^ 

(-1)^ K-1 - dT £^ + ( S 1 - E - ^^^+^+1) ( ) = 



k-2 



i=0 



(28) 
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where the terms in brackets ( ) contain relations involving partial derivatives of the La- 

grangian function L and applications of the total derivative dx, which are not written for simplic- 
ity. These equations may or may not be compatible, and a sufficient condition for compatibility 
is the regularity of the Lagrangian density C. In particular, we have: 



Proposition 5 If C £ Q^{J^-k) is a regular Lagrangian density, then there exists a unique 
vector field Xq € Xwi(Wo) which is a solution to equation (|2ip .- it is tangent to Wi, and is a 
semispray of type 1 in Wo- 



(Proof) As the Lagrangian density C is regular, the Hessian matrix — ^ — j is regular at 

every point, and this enables us to solve the above k systems of n equations (|28|) determining 
all the functions uniquely, as follows 

F,^ = qf+i , {k^i^2k-2) (29) 

In this way, the tangency condition holds for Xq at every point on Wi. Furthermore, the 
equalities show that Xq is a semispray of type 1 in Wo with local expression 

9 Ad d dt d , , s 9 

^° = « + a? + ^--'94:; ^WtW,^ ''^"M ■ 



However, if C is not regular, the equations ()28p may or may not be compatible, and the 
compatibility condition may give rise to new constraints. In the most favourable cases, there is a 
submanifold W/ Wi (it could be W/ = Wi) such that there exist vector fields Xo S Xw^ (Wo), 
tangent to W/, which are solutions to the equations 

i(Xo)Oo|w, =0 , 2(Xo)(pS)*??U, = 1 . (31) 
Finally, the relation among the results obtained in the two last sections is as follows: 



Theorem 1 The following assertions on a holonomic section ipo G r(yoj^) are equivalent: 

1. ipo is a solution to equation (jlip . that is, 

ro i{y)^o = 0, for every Y G X(Wo) . 

2. If Ipo is given locally by ipoit) = {t, qfit), qf{t),p\{t)), ^ « ^ /c - 1, /c ^ j ^ 2A; - 1, then 
the components of ipo satisfy equations (fH|) and (fT5|) . that is, the following system of kn 
differential equations 

^ <^^) 

3. ipo is a solution to the equation 

i(V';)(f^ooVo) = 0, (33) 
where ip'^: ^ TWo is the canonical lifting of V'o to the tangent bundle. 
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4- ipo is an integral curve of a vector field contained in a class of p^-transverse semisprays of 
type 1, {Xo} C X{yVo), satisfying the first equation in (fT9]) . that is, 

i{Xo)no = o. 



(Proof) 

(1 2) As we have seen in Section [3.2. H equation pT]) gives, in natural coordinates, the 
equations p^ . (fT5]) . ([T6]) and (fT7|) . As stated there, equation (fT3|) is redundant, since it 

is a combination of the others, and from equations (fT6]l we deduce that the section ipo G r(pj) 
hes in the submanifold Wi. Hence, equation pT]) is locally equivalent to equations (fT5]) . 
and (fT7|l . However, as we assume that is holonomic, equations (fT7|) hold identically, and thus 
equation ([TT]) is locally equivalent to equations and p!5|) . that is, to equations ([32]) . 



(2 <^ 3) If Vo(i) = {t, qt{t),qf{t),p\{t)) is the local expression of ijjo in natural coordinates, 
then V'ol*) = (li '7/(0iPa(^)); ^^i"^ t^i^ inner product i{il)'^){Vlo o gives, in coordinates, 

+ (^-^ - Pa ' - p'aJ dqf + (^p'X^^ - dq^ + {qf - qf+^) dp\ • 

Now, requiring this last expression to vanish, we obtain the system of {2k + l)n + 1 equations 

i -A -A • 4 „ .n dL ■ dL •_i 

PaC+i-C^+PaC+i = ; PA=g^ ; Pa = ^-Pa 

fc— 1 -A A 

Pa = ; 1i = Qi+i 



Observe that this system of equations is the same given by P^ , P^ , P^ , and p!7|) . 
The same remarks given in the proof of (1 44> 2) apply in this case. In particular, the fifth group 
of kn equations qf = q^^ is identically satisfied by the section since we assume it to be 
holonomic. Thus, bearing in mind the above item, we have proved that equation (|33p is locally 
equivalent to the kn differential equations (j32|) . 

(2 44> 4) As we have seen in this Section, if a generic vector field Xq S X(Wo) is given locally 
by p^ . then the first equation in ([19]) is locally equivalent to equations (f22]) . (f23]) . ([2^ and 
(|26p . As already stated, equation (j22p is redundant, since it is a combination of the others; and 
the n equations (|26p state, in coordinates, the result given in Proposition [3l In addition, since 
the vector fields Xq in the class are semisprays of type 1, the kn equations (|23p are identically 
satisfied. Thus, the first equation in pop is locally equivalent to the kn equations (|24p . Finally, 
the pj^-transverse condition for the class {Xo} is locally equivalent to / 7^ 0. 

Now, let cr G r(pj^) be an integral curve of Xq, that is, a' = Xq o a. If ex is given locally 
by a{t) = {t,qf{t),qf{t),p\{t)), then a'{t) = {l,qf{t),qf{t),p\{t)), and, taking / = 1 as a 
representative of the class {^ojj the condition of a to be an integral curve is locally equivalent 
to the equations 

qf = f^^°cT ; qf = Ffoa ; p\ = G\oa. 
Replacing these equations in (j23]) and (l24]l . we obtain the following 2kn differential equations 

1. =..+1 ; p'A=g^ ; PA = g^-PA ■ 
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Observe that, as every vector field in the class is 
are identically satisfied. Thus, the condition of 
semispray of type 1, Xq G XiWo), satisfying the 
equations ([32]) . 



a semispray of type 1, the first kn equations 
a to be an integral curve of a pj^-transverse 
first equation in (fT9l) is locally equivalent to 



4 Lagrangian formalism 
4.1 General setting 

Now we recover the Lagrangian dynamics from the unified formalism. We do not distinguish 
between the regular and singular cases, since the results remain the same in either case, but a 
few comments on the singular case will be given. First, we have: 

Proposition 6 The map Pi = p" ° ji ■ Wi — )• J'^^~^tt is a dijfeomorphism. 

(Proof) As Wi = graph. J^C, we have that J'^^^^tt ~ Wi. Furthermore, p\ is a surjective 
submersion and, by the equality between dimensions, it is also an injective immersion and hence 
it is a diffeomorphism. ■ 

Now, we must define the Poincare-Cartan forms in order to establish the dynamical equations 
for the Lagrangian formalism. First, we have the following result: 

Lemma 2 Let @k~i e {T* {J^-^ t:)) , Vl^^i = -d@k~i e Q'^{T*{J^-'^'k)) he the canonical 
forms in T*{J^^^tt). We define the Poincare-Cartan forms as Qc = ©fc-i G ^?^( J^'^'^vr), 
f]^ = -dQc = TC*nk-i G I?2(j2fc-i^)_ y/jg^ 0^ ^ {pi)*ec and Qo = {piT^C- 

(Proof) We have for Qc- 

{plTQc = [plTiTcQk-i) = {tI o piyQk-i = [pD^Qk-i = Qo , 

and for il^: 

{pl)*nc = {pl)*{-dec) = -d(p?)*0£ = -dOo = no • 



Alternatively, according to [12] and [l3] (see also [T], [26]), we can define the Poincare-Cartan 
1-form using the canonical structures of the higher-order jet bundles; in particular, 

Oc = Sj^^'^dL) +Ce I?^(J2^-V) , 

(k) 

where Sr) is the generalization to higher-order jet bundles of the operator used in the clas- 
sical Hamilton-Cartan formalism for problems in the calculus of variations which involve time 
explicitly (see [42] and [43] for details). 

Using natural coordinates, the local expression of the Poincare-Cartan 1-form is 

®^ = E E(-1)'^T ^ (dqti - q^dt) + Ldt . (34) 

r=li=0 \"^r+ij 



P.D. Prieto-Martmez, N. Roman-Roy: Unified formalism for higher-order non-autonomous. . .20 



Remark: 0£ is a 7r^^-|^^-semibasic 1-form. 

Prom the Poincare-Cartan 1-form, the concept of regularity for a higher-order Lagrangian 
density is a straightforward generahzation of the well-known definition for first-order non- 
autonomous dynamical systems. In fact, first we define the Poincare-Cartan 2-form as ^Ic = 
-dQc G ^?2(J2'=-^7^). Then 

Definition 10 A Lagrangian density C G f2^{J^TT) is regular if the Poincare-Cartan 2-form Qc 
has maximal rank. Elsewhere C is singular. 

In natural coordinates, the local expression of the 2-form is 

From this expression in local coordinates, we can see that the regularity condition for C is 
equivalent to 



for every y G J^'^^^vr. Thus, this notion of regularity is equivalent to the one given before. 
Geometrically, C is regular if, and only if, (Jl^, (7f^'^~^)*r/) is a cosymplectic structure on J^'^^^vr, 
that is, and (vf^^^"'^)*?] are both closed and Jl^" A (7f^^~^)*r/ is a volume form. 



4.2 Dynamical equations for sections 

Using the previous results, we can recover the Lagrangian sections in J^^~^7r from the sections 
in the unified formalism. 



Proposition 7 Let ipo £ r(pj^) be a holonomic section solution to equation (jlip . Then the 
section ipc = Pi ° 4'o ^ T{Tt'^^^^) is holonomic, and is a solution to the equation 

i^l i{Y)nc = 0, for every Y G X( J^'^" V) (36) 

{Proof) Since, by definition, ipo G ^(^1^) is holonomic if pi o Tp^ ^ r(7f^^~^) is holonomic, it is 
obvious that tl^c = Pi o V'o is a holonomic section. 

Now, recall that, since pi is a submersion, for every Y G X{J^^-'^tt) there exist some Z G 
X{yVo) such that pl^Z = Y, that is, Y and Z are pf-related. Note that this vector field is not 
unique, since Z -\- Z^, with Zq G ker/j°^ is also pf-related with Y. Thus, using this particular 
choice of p°-related vector fields, we have 

rc iiy)^c = {Pi o i^oT i{Y)nc = roMT i{Y)nc) = ro{i{z){pirnc) = roi{z)^o . 

Since the equality tpli{Z)Qo = holds for every Z G X{Wo), in particular it holds for every 
Z G X(Wo) which is p°-related with Y G X(J^^~^7r). Hence, we obtain 



rci{Y)Qc = roi{Z)no = 0. 
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The diagram for this situation is the following: 




Remark: Observe that, from this result, we have no equivalence between section ipo € r( 
solutions to equation ([TT]) and section il^c G r(7f^'^~^) solutions to equation ([36l) . but only that 
every holonomic section tpo solution to the dynamical equations in the unified formalism can be 
projected to a holonomic section ipc solution to the Lagrangian equations. Nevertheless, recall 
that section ipo solutions to equation take their values in the submanifold Wi, which is 
diffeomorphic to J^'^^^vr, and thus it is possible to establish an equivalence using the diffeomor- 
phism p\. 

Assume e r(/3°) is given locally by tpoit) = {t,qf{t),qf{t),p\{t)), < i ^ /c - 1, 
k ^ j ^ 2k — 1. Since ipo is assumed to be a holonomic section solution to equation ([TT]) . it 
must satisfy equations , ([T5|) and (fT7|) . The last group of equations is automatically satisfied 
because of the holonomy condition. Now, bearing in mind that the section ipo takes values in the 
submanifold Wi, and the characterization of Wi given in Proposition [21 equations ()14p and (jl5p 
can be /of-projected to J^'^~^7r, thus giving the following equations for the section '4>c = Pi° '^o- 



dL 



dqn 



d dL 



dt dq^ 



+ 



d2 dL 



dt^ dq^ 



+ ... + (-1) 



d'^ dL 



dt^ dq^ 



0. 



Finally, bearing in mind that '(/'£ is holonomic in J"^^' 
local expression is = {t,qQ{t)), such that j'^^~^4 
be rewritten in the following form 



TT, there exists a section (p G r(7r), whose 
= ipCj fiiid thus the above equations can 



dL 



dq^ 



d dL 



dt dqf 



+ 



d^ dL 



dt2 dq^ 



+ ... + {-!) 



d^ dL 



dt^ dqf 



0. 



(37) 



Therefore, we obtain the Euler-Lagrange equations for a kih. order non-autonomous system. As 
stated before, equation (j37p may or may not be compatible, and in this last case a constraint 
algorithm must be used in order to obtain a submanifold Sf ^ J'^^^^tt (if such submanifold 
exists) where the equations can be solved. 



4.3 Dynamical equations for vector fields 

Now, using the results stated at the beginning of the Section, we can recover a vector field 
solution to the Lagrangian equations starting from a vector field solution to the equation in the 
unified formalism. First we have: 



Lemma 3 Let Xq G X(Wo) he a vector field tangent to Wi. Then there exists a unique vector 
field Xc G je(j2^-i7r) such that o p° o = Tp" o o j^. 
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(Proof) Since Xq is tangent to Wi, there exists a vector field Xi G X(Wi) such that Tji oXi = 
° ji- Furthermore, as p\ is a diffeomorphism, there is a unique vector field Xc € X( J^'^^^vr) 
which is /o}-related with Xi; that is, Xc ° Pi = Tp} o Xi. Then 

Xcop^o ji = Xcop\ = Tp\ o Xi = Tpl o Tji o Xi = Tpl o Xo o ji . 



The above result states that for every Xo G Xwi(VVa) there exists a vector field X/; G 
X( J^'^~^7r) such that the following diagram commutes 



T(J^'=-%) 



J 



2A;-1 




As a consequence we obtain: 

Theorem 2 Lei Xq £ jCvv^ (VVo) a vector field solution to equations (|2ip and tangent to Wi 
f'ai /easi on i/ie points of a submanifold W/ M> Wij. T/ien i/iere exisis a unique semispray of 
type k, Xc G X( J^'^^-'^tt), which is a solution to the equations 



i{Xc)nc = ^ , iiXc)i7t^'-'rr^ = 1 



(38) 



(at least on the points of Sj = Pi(W/)/ In addition, if C € Q^{J^tt) is a regular Lagrangian 
density, then Xc is a semispray of type 1. 

Conversely, if Xc G X(J^^~^7r) is a semispray of type k (resp., of type I), which is a solution 
to equations (j38p (at least on the points of a submanifold Sf ^ J'^^~^'k), then there exists a 
unique vector field Xq G X>Vi(Wo) which is a solution to equations ()2ip (at least on the points 
ofWf = {p\)^^{Sf) ^ Wi ^ Woj, and it is a semispray of type k in Wo (resp., of type 1). 



{Proof) Applying Lemmas [2] and [3l we have: 

^=l{Xo)^o\^^ = i{Xo){pir^c\y^^ = {pir l{Xc)^c\yv, 



{pD* i{Xc){7:' 



1 = i{Xo){fJiYv\y,, = i{Xo){7r''~' o pi)*^ 
However, as p" is a surjective submersion, this is equivalent to 

= i{Xc)^c\pO(yv^) = ii.Xc)^c\j2k-i„ , 
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since /5i(Wi) = J'^^~^tt (or the submanifold Sj ^ J^'^^^vr). The converse is immediate, reversing 
this reasoning. 

In order to prove that Xc is a semispray of type k, we compute its local expression in 
coordinates. From the local expression (|27|) for the vector field Xq (where the functions Fj^ are 
the solutions of equations ([28|) ). and using Lemma O we obtain that the local expression of the 
vector field Xc G X{J^''-^7r) is 



1=0 * j=k J 

which is the local expression for a semispray of type k in J^^^^vr. 

Finally, if £ G Q^{J^tt) is a regular Lagrangian density, equations ([28|) become (p9|) . and 
hence the local expression of Xc is 



which is the local expression for a semispray of type 1 in J'^^~^'it. 
Remarks: 



• It is important to point out that, if £ is not a regular Lagrangian density, then Xq is a 
semispray of type k in Wo, but not necessarily a semispray of type 1. This means that 
Xc may be a solution to the Lagrangian equations for vector fields, but the trajectories 
given by its integral sections are not solutions to the dynamical system (the sections 
solution to the dynamical problem must be holonomic, but the integral sections of Xc 
are only holonomic of type k). Thus, for singular Lagrangians, this must be imposed 
as an additional condition. This constitutes a relevant difference from the case of first- 
order dynamical systems, where this condition {Xc is a semispray of type 1) is obtained 
straightforwardly in the unified formalism. 

For singular Lagrangians, only in the most interesting cases can we assure the existence 
of a submanifold W/ ^ Wi and vector fields Xq G (Wo) tangent to W/ which are 
solutions to equations ([3T]) . Then, considering the submanifold Sf = p\iyVf) ^ J^'^~^7r, 
in the best cases we have that those semisprays of type 1 Xc exist, perhaps on another 
submanifold ^ Sf where they are tangent, and are solutions to equations 

i{Xc)nc\M=^ , z(X^)(7f2^-i)*r? =1. (39) 

^ M f 

• Notice that Theorem [2] states that there is a one-to-one correspondence between vector 
field Xo G X>Vj(Wo) solutions to equations ([2T|) and vector field Xc G X{J'^^~^7t) solutions 
to ([38l) . but not uniqueness. In fact, we cannot assure uniqueness of the vector field Xc 
unless the Lagrangian density is regular, as we can see in the following result: 

Corollary 2 If the Lagrangian density C G Q^{J^7r) is regular, then there is a unique semispray 
of type I, Xc G X( J^'^^-'^vr), which is a solution to equations (pHj) . 
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(Proof) If the Lagrangian density C £ Q^{.J^ vr) is regular, using Proposition O there exists a 
unique semispray of type 1, Xq S j£(Wo), solution to equations (f2T]) and tangent to Wi. Then, 
using Theorem [21 there is a unique vector field Xc G X(j2^-V), which is a semispray of type 1 
in J^'^~^7r and is a solution to equations (p8|) . ■ 

In other words, uniqueness of the vector field Xc is a consequence of uniqueness of Xq. 
Finally, as a consequence of Theorem [1] and the results stated in this Section, we obtain: 

Theorem 3 The following assertions on a section (j) G r(7r) are equivalent: 

1- CL solution to equation (j36p . that is, 

(fk-i^y = 0, for every Y G X{J^''-^tt) . 



2. In natural coordinates, if (f) = (t, go(t)), then j (j) = {'t,Qoi^)'Qii^)^---^l2k-i(^)) ^■^ ^ 
solution to the kth order Euler- Lagrange equations given by ([ST]) , that is. 



dL 



d dL 



dt dqf 



+ 



d2 dL 



dt2 dq^ 



+ ... + {-!) 



k d'' dL 



dt^' dq^ 



0. 



3. Denoting ipc = j"^^ ^<P, then il^c is a solution to the equation 

i{^p'c){nco4^c) = o, 

where M — )• T{J'^^~^tt) is the canonical lifting of ijjc to the tangent bundle. 



4. ^4> is an integral curve of a vector field contained in a class of t: 



2k~l 



transverse semis- 



prays of type \, {Xc} C X(J vr), satisfying the first equation in (pSl) . that is. 



l{Xc)^L 



0. 



5 Hamiltonian formalism 

5.1 General setting 

In order to describe the Hamiltonian formalism on the basis of the unified one, we must distin- 
guish between the regular and non-regular cases. In fact, the only "non-regular" case we consider 
is the almost-regular one, so we need to define the concept of almost-regular Lagrangian density. 

Before doing so, we must define the generalization of the Legendre map from the first-order 
time-dependent case. Since @c G i7^( J^'^^^vr) is a 7r^'j|~[^-semibasic 1-form, we can give the 
following definition: 



Definition 11 The extended Legendre-Ostrogradsky map associated with the Lagrangian den- 
sity C is the map TL: j'^^^^ix — t- T*(J^~^7r) defined as follows: for every u G T( j2'^'~"^7r), 

G^n) = (T^f_-^n) | tI{tj.u-^M)) , 

where Tj2k~i^: T{J'^^^'^tt) — J'^^^^-k is the canonical submersion. 
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This map verifies that jk-i^oT L = t^I^Li^, where TTjk-i^ : T*{J^ ^vr) — ^vr is the natural 
projection. Furthermore, if efc_i G I?i(T*( J^^-^Tr)) and Qj^^i = -de^-i G n'^{T*{J^-^7r)) are 
the canonical 1 and 2 forms of the cotangent bundle T*(j'^^^7r), we have that 

J^C Qk-l = @c ! ^k-i = ■ 



Bearing in mind the local expression @ of the tautological 1-form on T*{J^ ^tt) and the 
local expression (j34p of 0£, we have that the local expression of the map J^C is: 



that is, this map coincides with the extended Legendre-Ostrogradsky map defined locally in 
Section [3.2.H thus justifying the notation and terminology introduced therein. 

Notice that dimT*(j'=~i7r) = 2A;n + 2 > 2A;n + 1 = dim J^'^-^vr. Thus, T*(j''-^7r) is not a 
suitable dual bundle to J'^^~^tt for giving a Hamiltonian description of the dynamical system. 
Therefore, according to, for instance, [41] and the references therein, we consider the bundle 
jk-i^* ^ T*(j'=-i7r)/(7f'=-i)*T*R, with the natural projections 

rp*/ jk—l„\ . jfc— 1„* „r . jk—l„* . jk—l„ - r „ =.k—l . jk—l„* . m> 

fi: L [J vrj J vr , tTj^-i^: J vr J vr , r = i^jk-i^ o tt : J vr M, 

where 7r^fe-i^ is the map satisfying -Kjk-i^ = TTjk-i^ ° fJ-. Notice that dim j'^^-'^vr* = 2kn + 1. 

Thus, we define the restricted Legendre-Ostrogradsky map as TC = fi o TL : J^'^^^tt — )• 
jk-i^* ^ This map satisfies Tr^j..^^ o J^C = tt^^^^, and has the following local expression 



J'CH = t 



k—r / 
i=0 \ 



dL 

dqf^ 



r+i 



In other words, this map coincides with the restricted Legendre-Ostrogradsky map defined locally 
in Section [3.2.11 This justifies the notation and terminology introduced in that Section. 



Proposition 8 For every y G J^'^ ^tt we have that rank(J^£(y)) = rank(J^£(?/)). 



We do not prove this result. Following the patterns in [TT], the idea is to compute in natural 
coordinates the local expressions of the Jacobian matrices of J-'C and J-C Then, observe that 
the ranks of both maps depend on the rank of the Hessian matrix of L with respect to q^ at 
the point y, and that the additional row in the Jacobian matrix of J-C is a linear combination 
of the others. See [17] for details in the first-order case. 

As a consequence of Proposition [8l and taking into account the different definitions given for 
the regularity of the Lagrangian density, we arrive at the following result: 

Proposition 9 Given a Lagrangian C G f2^{J^'K), the following statements are equivalent: 

1. Vtc has maximal rank on j'^^~^Tr. 

2. The pair (ri^, (yf^'^^^)*??) is a cosymplectic structure on J'^'^^^n. 
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3. J-C: J"^^ "^TT —7- ^vr* is a local diffeomorphism. 

4. TL: J'^'^^^TT — )• T*(j'^~^7r) is an immersion. 



(Proof) It is easy to check that ah the statements are locahy equivalent to 

Mkdqi 



Now, we denote hy V = lm{FC) = TC{J^^-^^) A T*{J^-^^) the image of the extended 

Legendre-Ostrogradsky map; and hyV = lm{TC) = TC{j'^^-^n) 4 J^-^tt* the image of the 
restricted Legendre-Ostrogradsky map. Let Tq = f o j: V ^ ^ he the natural projection. We 
can now give the following definition: 

Definition 12 A Lagrangian C G Q^{J^'k) is called an almost-regular Lagrangian density if: 

1. V is a closed submanifold of J^^'^tt* . 

2. J-C is a submersion onto its image. 

3. For every y G J'^''~^tt, the fibers TC'~^{TC{y)) are connected submanifolds 0/ J^'^~^7r. 

As a consequence of Prop. [HI we have that V is diffeomorphic to V. This diffeomorphism is 
just /i restricted to the image set V, and we denote it by Jl. This enables us to state: 

Lemma 4 If the Lagrangian density C G fi^{J^'K) is, at least, almost-regular, the Hamiltonian 
section h G r(/i>v) induces a Hamiltonian section h G r(/i) defined by 

h{[a]) = (p2 o hmpl)-\j{[a]m, for every [a] G V. (40) 

{Proof) It is clear that, given [a] G J^~^Tr*, the section h maps every point (y, [a]) G (P2)~"^(['^]) 
into P2^[p2{h{y, [a]))]- So we have the diagram 



3 rp* / jk—1 




Thus, the crucial point is the /72-projectability of the local function H. However, since a local 
base for ker p2* is given by 

/ d d \ 

ker p2* 



we have that H is /02-pi'ojectable if and only if 

This condition is fulfilled when [a] G V, which implies that /02[^((P2) "'^(['^]))] ■ 
Remark: In the hyperregular case, we have V = J^^^ir*. 

Locally, this Hamiltonian |U-section is specified by the local Hamiltonian function H G 
C°°(J*^-i7r*), that is, 

h{t,qf,p\) = {t,qf,-H,p\). 
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5.2 Hyperregular and regular systems. Dynamical equations for sections and 
vector fields 

Now we analyze the case when £ is a regular Lagrangian density, although by simplicity we 
focus on the hyperregular case (the regular case is recovered from this by restriction on the 
corresponding open sets where J-^C is a local diffeomorphism) . This means that the phase space 
of the system is J^^^ir* (or the corresponding open sets). 

In this case, we can give the explicit expression for the local Hamiltonian function, which is 

H = Y^PAQti + p\-\FC-'Yqt - (vrf ~i o ^£~i)*L . (41) 
1=0 

The Hamiltonian section h is used to construct the Hamilton-Cartan forms in J^^^tt* by 
making 

eh = h*ek-i€ f2\j''-^7r*) , nh = h*nk^ie n^{j''-^7r*), 

where Qk-i and ^k-i are the canonical 1 and 2 forms of the cotangent bundle T*(j''~^7r). 
Bearing in mind the local expression of @k-i and f^fe-i, the local expression of the forms Qh 
and ri/i is 

= p\dqf - Hdt , Uh = dqf A dp\ + dH A dt , 
Notice that J='C*eh = ©£ and J='C*nh = ^c- 

Proposition 10 If C £ Q^{J'''k) is a hyperregular Lagrangian, then = P2°ji '■ VVi — J'^ ""^vr* 
is a diffeomorphism. 



(Proof) The following diagram is commutative 




that is, we have = P2 ° ji = J^Co p\. Now, by Proposition [6l the map p\ is a diffeomorphism. 
In addition, as C is hyperregular, the map J^C is also a diffeomorphism, and thus P2 is a 
composition of diffeomorphisms, and hence a diffeomorphism itself. ■ 

This last result allows us to recover the Hamiltonian formalism in the same way we recovered 
the Lagrangian one (see Section H]), just using the diffeomorphism to define a correspondence 
between the solutions of both equations. 

Using the previous results, we can recover the Hamiltonian sections in J^^^tt* from the 
sections solution to the equations in the unified formalism. 

Proposition 11 Let C G Q^{J^tt) he a hyperregular Lagrangian. Let ipo G ^(Pm) section 
solution to equation (llip . Then the section = P2 ° i^o ^ ^(f) is a solution to the equation 



il^l i{Y)Vth = 0, for every Y G X( j''" V*) 



(42) 
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(Proof) The proof of this result is analogous to the proof given for Proposition d 
The diagram for this situation is the following: 

Wo. 




Remarks: 



• Observe that, for the Hamiltonian sections, the condition of holonomy on the section V'o 
is not required. This is because we only need tpo to be a holonomic section of type /c, and 
this condition is always fulfilled. 

• As for the Lagrangian sections given by Proposition [71 this last result does not give an 
equivalence between sections V'o £ r(p^), which are solutions to equation (fTT]l . and sections 
iph £ r(f ), which are solutions to equation ()42p . However, recall that sections which 
are solutions to the dynamical equations in the unified formalism, take values in Wi, and 
hence we are able to establish the equivalence using the diffeomorphism 



Let Vo(i) = {t,qf{t),qf{t),p\{t)) G T{p'^), ^ i ^ A; - 1, A; < j ^ 2A; - 1, be a solution to 
equation (fTTl) . Hence, iIjq must satisfy equations (fT4]l . (fT5]) and (fTT]) . Now, bearing in mind the 
local expression for the local Hamiltonian function H given in (|4ip . we obtain the following 2kn 
equations for the section iph = P2 ° i'o = {t,Qf{t),p\{t)): 



dH 



dH 



Pa 



dqi 



(43) 



So we obtain the Hamilton equations for a /cth-order non-autonomous system. 



Next, we recover the Hamiltonian vector field from the vector field solution to the dynamical 
equations (jl9p in the hyperregular case. As P2 is a diffeomorphism by Proposition [TOl the 
reasoning we follow is the same as that for the Lagrangian formalism. 



Lemma 5 Let C G Q^{J^'k) he a hyperregular Lagrangian. Let Xq G X(VVo) he a vector field 
tangent to Wi. Then there exists a unique vector field G X{J^~^tt*) such that X^ o P2° ji = 
Tpl o XoO ji . 

(Proof) The proof of this result is similar to the proof given for Lemma O ■ 
This result states that, for every Xq G 3£>Vi(Wo), we have a vector field X^ G X{j'^^^7r*) 
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such that the following diagram commutes 




Xh 



Theorem 4 Let C G f^^^j'^ir) be a hyperregular Lagrangian, and Xq G ^Wi(Wo) th^ vector 
field solution to equations ()2ip and tangent to Wi. Then, there exists a unique vector field 
Xh G 3i{J^~^'K*), which is a solution to the equations 



i{Xh)^h = ^ , i{Xh)f*7J = l. 



(44) 



Conversely, if Xh G -^vr*) is a solution to equations (I44p . then there exists a unique 

vector field Xq G X>v^(Wo), tangent to Wi, which is a solution to equations (|2ip . 



(Proof) The proof of this result is analogous to the first part of the proof given for Theorem [21 
Lemma [5] now being used to obtain the vector field Xh G X{J^~^7r*). ■ 

In local coordinates, if the vector field Xo G Xw^ (Wo) solution to equations (f2T]) is given by 
(j30p . by using Lemma [5] we obtain the local expression for the vector field Xh, which is 



d 



d dL d 

+ 



+ drip: 



'A) 



d 

dpj 



Finally, to close the hyperregular case, as a consequence of Theorem [T] and the results stated 
in this Section, we obtain the following result: 



Theorem 5 The following assertions on a section iph £ '^{t) fi^e equivalent: 



1. iph is a solution to equation (j42|) . that is, 

i{Y)nh = 0, for every Y G X( J^" V*) . 

2. In natural coordinates, if iph is given by iPhit) = {t,qf{t),p\{t)), ^ i ^ - 1, then the 
components of iph satisfy the kth order Hamilton equations given by ()43p . that is. 



dH 

5p* 



p'a 



■>Ph 



dH 
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3. iph is a solution to the equation 

i{'^'h){^h o ^h) = 0, 
where ip'f^:M.^ T{J^^^7r*) is the canonical lifting of tph to the tangent bundle. 

4- iph is an integral curve of a vector field contained in a class of f-transverse vector fields, 
{Xh} C X(j'^~^7r*), satisfying the first equation in (HH), that is, 

iiXh)nh = o. 

5.3 Singular (almost-regular) Lagrangians. Dynamical equations for sections 
and vector fields 

Recall that, for almost-regular Lagrangians, only in the most favourable cases can we assure the 
existence of some submanifold W/ ^ Wi where the dynamical equations can be solved. In this 
case, the solutions to the Hamiltonian formalism cannot be obtained straightforwardly from the 
solutions in the unified formalism, but rather by passing through the Lagrangian formalism and 
using the Legendre-Ostrogradsky map. 

In this case, the phase space of the system is V = lm{TC) J'^^^tt*. We denote by 
TCo'. J^'^^^vr — V the map defined by = j o TLq- As in the hyperregular case, the 
Hamiltonian section h is used to construct the Hamilton-Cartan forms on V as follows: 

They verify that -F£*e^ = Oc and TC^Ql = Qc- 

Proposition 12 Let C G Q^{J'^it) he an almost-regular Lagrangian. Letxpo G r(/9g) be a section 
solution to equation (llip . Then, the section ip'^ = TCoOipc £ r(ro) is a solution to the equation 

«)*i(y)f]^ = 0, for every Y e?i{V). (45) 

(Proof) Since the Lagrangian density is almost-regular, the map J- Co is a submersion onto 
its image, V. Hence, for every Y G X{V) there exist some Z G X^J'^^^^tt) such that Z is 
-F£o-related with Y, that is, TCq^Z = Y. Using this, we have 

Then, using Proposition [71 we have proved 

{rhri{YWh = rci{mc = ^. 

■ 

The diagram for this situation is the following: 
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Now, assume that there exists a submanifold W/ ^ Wi and vector fields Xq G (Wo) 
tangent to W/ which are solutions to equations ([3T]) . Now consider the submanifolds 5"/ = 
^ J^'^-^TT and P/ = pliWf) = FC{Sf) ^ J''~^tt*. Using Theorem El from the 

vector fields Xq € ^Wi(Wo) we obtain the corresponding vector fields Xc G X{J'^^~^7r), and 
from these, the semisprays of type 1 (if they exist), which are perhaps defined on a submanifold 
Mf ^ Sf, are tangent to Adj and are solutions to equations (j39p . So we have the diagram 




Now, following analogous procedures for autonomous and non-autonomous systems jl8l I29j . 
one can prove that there are semisprays of type 1 in Mf (perhaps only on the points of another 
submanifold Mf ^ Mf), which are J^£-projectable on Pf. These vector fields X^ = TL^Xc G 
j£(P) are tangent to Pf and are solutions to equations 

i{Xl)^l\p^=^ , i[Xl)flTt\p^ = \ (46) 

Conversely, as TCq is a submersion, for every vector field G X(7^) solution to equations 
dMD, there is a semispray of type 1, X^ G X(J^'^~-^7r), such that FCq^Xc = Xf^, and we can 
recover solutions to equations (f3T]) using Theorem [2l 

Of course, for the almost-regular case, we have a similar result to Theorem [5l on the points 
of the final constraint submanifold Pf. 

6 Examples 

6.1 The shape of a deformed elastic cyhndrical beam with fixed ends 

As a first example we consider a deformed elastic cylindrical beam with both ends fixed. The 
problem is to determinate its shape; that is, the width of every section transversal to the axis. 
This system has been studied on many occasions, such as \9] (Chapter 3, §3.9) and [25] (Chapter 
IV, §4). Strictly speaking, it is not a time-dependent mechanical system, but it can be modeled 
using a configuration bundle over a compact subset of M, where the base coordinate represents 
every transversal section of the beam, thus allowing us to show an application of our formalism. 
For simplicity, instead of a compact subset, we take the whole real line as the base manifold. 
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The configuration bundle for this system is vr: E ^ where E is a 2-dimensional smooth 
manifold. Let x be the global coordinate in M, and t] G i7^(R) the volume form in M with local 
expression rj = dx. Natural coordinates in E adapted to the bundle structure are {x,qo). Now, 
taking natural coordinates in the higher-order jet bundle of vr, the second-order Lagrangian 
density for this system, C G i7^(J^7r), is locally given by 

£{x,qo,qi,q2) = L- {Tt'^)*ri= {]^p{x)ql + p{x)q^ dx , 

where p, p ^ C°°( J^vr) are functions that only depend on the coordinate x and represent physical 
parameters of the beam: p is the linear density and p is a non-vanishing function involving 
Young's modulus of the material, the radius of curvature and the sectional moment of the cross- 
section considered (see [9] for a detailed description) . This is a regular Lagrangian density, since 
the Hessian matrix of the Lagrangian function L G C°°(J^7r) associated with £ with respect to 
92 is 

and this 1x1 matrix has maximum rank, since p is a non-vanishing function. 

Remark: If the beam is homogeneous, p and p are constants (with p ^ 0), and thus the 
Lagrangian density is "autonomous" , that is, it does not depend explicitly on the coordinate of 
the base manifold. This case is analyzed in [25j . 

As this is a second-order system, we consider the bundles W = J^n x ji^ T*{J^Tr) and Wr = 
J^ir Xji^ J^TT*, with natural coordinates {x,qo,qi,q2,q3,p,p^,p^) and ix,qo,qi,q2,q3,p'^,P^), 
respectively. Now, using the notation and terminology introduced throughout this article, if 
Qi G f?i(T*(J^7r)) and Qi G i7'^{T*{J^TT)) are the canonical forms of T*{J^tt), we define the 
forms = P2&1 G 17^(VV) and Q = P2^i G i7^(W), whose local expressions are 

Q = p'^dqo + p^dqi + pdx ; Q = dqo A dp^ + dqi A dp^ - dp A dx . 
The coupling 1-form C G i7^(W) has the local expression 

C = C ■ p^r] = {p + p^qi + p^q2)dx , 
and then we can introduce the Hamiltonian submanifold 

Wo = [we W: £{w) =C{w)^ 4 W, 
which is locally defined by the constraint function C — L = 0, whose coordinate expression is 

C - L=p + p^qi + p^q2 - ^p{x)ql - p{x)qo = . 

Finally, we construct the Hamiltonian /iy\;-section h G T{py^;), which is specified by giving the 
local Hamiltonian function H, whose local expression is 

H{x,qo,qi,q2,q3,P^,P^) = P^qi +P^q2 - ^p{x)q2 - p{x)qo ; 

that is, we have h{x,qo,qi,q2,q3,p^ ,p^) = {x,qQ,qi,q2,q3,—H,p'^,p^). Using this Hamiltonian 
section, we define the forms Go = j*<d G f2^{yVo) and Oo = j*0 G ^^^(Wo), with local expressions 

Go = p^dqo + p^dqi + Q^(x)g| + p{x)qo - p^qi - p^q2^ dx , 

= dqo A dp^ + dqi A dp^ + [-p{x)dqo + p^dqi + {p^ - p{x)q2)dq2 + qidp^ + q2dp^) A dx . 
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In order to state the Lagrangian-Hamiltonian problem for sections in this system, let Y G 
X(yVo) be a generic vector field locally given by 

ox oqo oqi oq2 oqs op op'- 

Now, \iipo{x) = {x,qo{x),qi{x),q2{x),qs{x),p^{x),p^{x)) is a holonomic section of the projection 
pj^, equation (fTTI) leads to the following 5 equations (the redundant equation (fTSl) is omitted): 

qo = qi ; qi = 52 (47) 

pO = p{x) ; p^ = (48) 
p^ = q2n{x) (49) 

Equations (j47p give us the condition of holonomy of type 2 for the section, which are also redun- 
dant since we assume that V'o is holonomic. Equation ()49p is a pointwise algebraic condition, 
from which we know that the section ij^o must lie in a submanifold Wi that can be identified 
with the graph of the extended Legendre-Ostrogradsky map, FC 

Now we compute the local expression of the map J-C: J^n — t- T*(J^7r). From Corollary [T] 
we know the general expression for this map, and we obtain: 

FC p° = -92^ - ; FC p^ = q2fi ; FC p = --^ml + ^1^2^ + gigs/" + qop ■ (50) 

Therefore, the section ipo G r(pj) is a holonomic section of the projection pj^, which lies in the 
submanifold Wi ^ Wo defined by the above constraint functions, and whose last components 
satisfy the differential equations 

p^ = p{x) ; p^ = —p^ . 

Now we state the Lagrangian-Hamiltonian problem for vector fields: we wish to find Xq G 
X(>Vo) solution to (fT9]) . If Xo is locally given by 

d d d d d d d 

Xo = 1- /07^ 1- /l7^ 1- F2-Z h F-^-- h G + G TT—j- , 

ox oqo oqi oq2 oqs op op 
then equations ([T9]) lead to the following (again, the redundant equation ([22]) is omitted): 

fo = f-qi ; fi = f-q2 (51) 

G^ = f-p{x) ■ G' = -f-p'' (52) 

/ = 1 (53) 

/ • (/ - q2^^{x)) = (54) 

Equations ([5T]) give us the condition of semispray of type 2 in Wo for Xo. In addition, equation 
()54p is an algebraic relation from which we obtain, in coordinates, the result stated in Proposi- 
tions [3] and HI that is, the vector field Xo is defined along a submanifold Wi which we identify 
with the graph of the extended Legendre-Ostrogradsky map and is defined by 

Wi = {weWo: Uw) = Uw) = 0} , 

where ir = P^ — FC p^ , r = 1,2. Thus, using ([5T]) . ([52|) and ([53]) . Xo is given locally by 



X — + + ^+f'^+F^+ (55) 

° dx dqa '^^ dqi ^ dq2 ^ dqs ^ dp^ ^ dp^ 
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Notice that the functions F2 and F3 in (j55|) are not determined until the tangency of Xq on Wi 
is required. This condition is locally equivalent to checking if the following identities hold 

UXo)Co\w,=0 , L(X„)6lw, =0. 

As we have seen in Section 13.2.21 these equations lead to the Lagrangian equations for the vector 
field Xq] that is, on the points of Wo we obtain 

L{Xo)^0 = P + + '^3^^ + + ^3/^ = ° 

L{Xo)^i = (gs - F2)fi = . (57) 

Equation (f57|) gives us the condition of semispray of type 1 for the vector field X^ (recall that 
is non- vanishing) , and equation (j56|) is the Euler-Lagrange equation for Xq. Observe that, 
since /i is a non- vanishing function, these equations have a unique solution for F2 and -F3. 
Hence, there is a unique vector field Xq G X(Wo) solution to the equations i{Xo)^o\yv-^ = and 
iiXo){p^)*rj\y^^ = 1, which is tangent to the submanifold Wi ^ Wq, and is given locally by 

„ 9 d d d 1 f d'^fi dfi\ d d r. d 
Xo = i:r- + qiT: hg27^ ^ Qst, P + Q2ir^ + 2^3 7^ 7^ ^ PTTn ~P 



ax ago oqi oq2 p- \ oxox ox J oq^ op" op^ 

Finally, we recover the Lagrangian and Hamiltonian solutions for sections and vector fields. 
For the Lagrangian solutions, by Proposition [71 from the holonomic section S ^(Pm) 
recover a holonomic section ipc = p° o ip^ ^ ^{^^) solution to equation p6p . In particular, if 
ipo{x) = {x,qo{x),qi{x),q2{x),q-i{x),p^{x),p^{x)), then ^pc{x) = {x,qo{x),qi{x),q2{x),qs{x)) is 
a holonomic section solution to equations (I48p . which, bearing in mind the local expression ()50|) 
of the extended Legendre-Ostrogradsky map, can be written locally as 

^ + ^2^ + ^2^-^+93/^ + 93^ =0 (58) 

ox Oxox ox 
(92 - 93)^ = (59) 

Equation ()59p gives the condition for the section ■0£ to be holonomic, and it is redundant since 
we required this condition to be fulfilled at the beginning. Now, if (j){x) = {x,y{x)) is a section 
of vr such that j^cp = if^c, then the Euler-Lagrange equation can be written locally 

-—{Hy)+p = {) . 
dx^ 

In the case of an homogeneous beam, the Euler-Lagrange equation reduces to /iy^™) + p = 

For the Lagrangian vector field, from Lemma [3] and Theorem [2l we can recover, from the 
semispray of type 1 Xq S X(Wo) a semispray of type 1, Xc G X(J^7r), which is a solution to 
equations ([38]) . and is locally given by 

d d d d I f d'^p. dp\ d 

= + ^17^ ^ ^27^ h 9375 [P + 927rir + ^^37^ -7^ ■ 

OX oqo oqi oq2 p \ oxox ox J oqs 

Now, as £ is a regular Lagrangian density, for the Hamiltonian solutions we can use the results 
stated in Section [5.21 and recover the Hamiltonian solutions directly from the unified formalism. 
For the Hamiltonian sections, using Proposition [Til from a section ipo £ r(/OK) fulfilling equation 
pT]) we can recover a section tph = P2 ° 'i'o ^ ^i^) solution to equation ([12]) . In particular, if 
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Tpoix) = {x,qo{x),qi{x),q2{x),q3{x),p^{x),p^{x)), then 7ph{x) = {x,qo{x),qi{x),p^{x),p^{x)) is 
a section solution to equations (|47p and (j48p . which can be written locally as 



^0 



dH 



dH 

dp^ 



p 



■>Ph 



dH 

dqo 



p 



dH 

dqi 



where H G C°°(J vr*) is the local Hamiltonian function with local expression 

H = p%^ + - pqo . 



2fi 



For the Hamiltonian vector field, from Lemma[S]and TheoremlU the vector field Xq G X(VVo) 
gives a vector field Xfi S X(J^7r*) solution to equations dM]), which is locally given by 



_ d d p^ d d 

^ dx dqo fi dqi ^ dp^ 



p 



_d_ 

dp^ 



6.2 The second-order relativistic particle subjected to a potential 

Consider a relativistic particle whose action is proportional to its extrinsic curvature [361 ESI 
[Ml Eh]. Now assume this system is subjected to the action of a generic potential depending on 
the time and the position of the particle, thus obtaining a time-dependent dynamical system. 

The configuration bundle for this system is M, where is a (n-|- l)-dimensional smooth 
manifold. Let t be the global coordinate in M, and r] G f?^(M) the volume form in M with local 
expression rj = dt. Natural coordinates in E adapted to the bundle structure are denoted by 
{t,qQ), 1 ^ i ^ n. Now, bearing in mind the natural coordinates in the higher-order jet bundle 
of vr, the second-order Lagrangian density for this system, £ G f2^{j'^Tr), is locally given by 

C{t, ql qU\) = Hf{q\f - {q\q\?] + V{t, g^)) dt = {j^^^J + ^(^^ ^l)) , 

where a is some nonzero constant and V G C°°(J^7r) is a function depending only on t and ^q. 
This is a singular Lagrangian density, as we can see by computing the Hessian matrix of the 
Lagrangian function L G C°°(J^7r) associated with C with respect to q2, which is 



d'^L 

dq^dqi 



—^mq\?-2{<f,)\<t,f)qfqt 

+{q{f{q{q\){q2qt - q?q2) - {ql?iqifq2q2] '^is^A 
-^[9- (qlfqiqt + mq2)qtq2 - iq\fq2q2] iiB = A, 



and a long calculation shows that det ( ^ . ) = 0. 



d^L 
dq^dq'2 

Consider the bundles W = J 27rXji^T*(jV) and Wr = vr X jijr J^vr*, with natural coordi- 
nates it,ql,q\,qi,qi,p,p'^i,p}) and it,qi),q\,qi,qi,p°i,pj), respectively Now, if Gi G f2^{T*{J^TT)) 
and ill G i7^(T*( J^vr)) are the canonical forms of the cotangent bundle of J^vr, we define 

e = p^ei = Pidql+p}dq\+pdt G I^^W) ; ft = pl^i = dqlAdp^i+dq{Adp}-dpAdt G n^{W) . 

The coupling 1-form C G f2^{W) has the local expression C = C ■ p^ij = {p + p^q\ +p}q2)dt, and 

from this we can introduce the Hamiltonian submanifold Wo '-^ W, which is locally defined by 
the constraint function C — L = 0, whose coordinate expression is 

C-L=p + p^qi + p}qi - J^Va - V{t, qi)) . 
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This allows us to construct the Hamiltonian ^yy-section h £ r{fj.y\;), which is specified by giving 
the local Hamiltonian function H, whose local expression is 



a 



V9-V(,t,i 



that is, we have h{t,qQ.,q\,q\,q\,p^,p\) = {t,qQ,q\,q2,ql,—H,p^,pl). Using this Hamiltonian 
section, we define the forms Qo = jo© ^ ^^O^o) and = j*n G /^^(Wo), with local expressions 



©o = P^dqi, + pldq{ V{t, qi,) - p^ " pIq^ 



= dql A dpy + dq\ A dpi + ( q^dp^ + Q#dp^ - — dg, 
a 



dV 



+ 



P\ 



dgf 



dq^ A dt . 



In order to state the Lagrangian-Hamiltonian problem for sections for this second-order 
system, let Y G X(Wo) be a generic vector field locally given by 



d 



.1 d 



aJ_ 

dt ' -"^ dq^ ' dqf ' ^ ' dq^ ' * dq^ ' " dp\ ' " dp\ 



Now, if ^o(i) = {t,qQ{t)-,q\{t)-,q2{t)-i'l\{t)-iP^{t)-,Pi{i)) is a holonomic section of the projection 
equation (fTT]l leads to the following 5n equations (the redundant equation (fT3]) is omitted): 



Qo=<lt 



P\ 



P\ 



dV 



P\ 



-P\ 



a 



a 



(60) 

[{{cAfW2f - 2{Mf) qt + {l<fi)Wifqt] (ei) 

(62) 



Equations ()60p give the condition of holonomy of type 2 for the section ipo, which are also redun- 
dant since the holonomy of tpo is already assumed. Equations (I62p are an algebraic condition, 
from which we conclude that the section ipo must lie in a submanifold Wi that can be identified 
with the graph of the extended Legendre-Ostrogradsky map, J-C. The expression in natural 
coordinates of this map FC: J^n — >• T*(J^7r) is obtained from Corollary [1] and is 



J'J^*p"a 
FC*p = 



a 



(9i)^vV 



+ 



a 



a 



imffiQkl) - iq\)\q\q2){qiqi) - (^i^)^) q^ - (gi) 



iq\rV9 



a 



iq\)W9^ 

M?q^ - {q\qi)qt 



(63) 



(^1 



a 



+ (q{fiqifiq{qi) - iq\?iq{qiMqi)) qt 



+- 



a 



{q\Y\n? 

a 



mq\f)\q\ql) - {<ii)\^Am<A) - {<Aqi)9) qt - WifaqtVi qt 



Wifqt-{q\q\)qt]qt 
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Hence, the section G '^(Pr) is holonomic and lies in the submanifold Wi ^ Wo defined by 
the constraint functions given above, and its last components satisfy the 2n differential equations 

p'a = ^ , p\ = -p'a- J(^r^ [m\lhf - it + iQ\Q2)iQ\)'4] ■ 

Now we state the Lagrangian-Hamiltonian problem for vector fields, that is, we wish to find 
a vector field Xq G X(VVo) solution to equations (fT9]) . If the vector field Xq is locally given by 

The equations ([T9]) lead to the following 5n + 1 equations (the redundant equation (j22|) is 
omitted) : 

fi = f-qt ; ft = f-<i2 (64) 

G^ = ^ ; G\ = -p\ - Wi?{<i^? - ^{<A<i2?) it + {^iq\M?4] (65) 

/ = 1 (66) 
f-{p\- Uf<l2 - WilD^t)) = (67) 



From equations (I64p we obtain the condition of semispray of type 2 for the vector field Xq. In 
addition, equations (j67p are algebraic relations between the coordinates in Wo which give, in 
coordinates, the result stated in Propositions [3] and HJ that is, the vector field Xq is defined along 
a submanifold Wi which we identify with the graph of the extended Legendre-Ostrogradsky. 
Thus, using (f64|l . (|65]) and (|66]l . the vector field is given locally by 



where the functions G\ are determined by (j65p . Since we wish to recover the solutions in the 
Lagrangian formalism from the vector field Xq, we must require it to be a semispray of type 1. 
This condition reduces the set of vector fields Xo G XiWo) given by ([68]) to the following ones 

V _ d A ^ A 9 A 9 TpA 9 dV d d , , 

Notice that the functions are not determinated until the tangency of the vector field Xq on 
Wi is required. From the expression in local coordinates (j63p of the map J-"£, we obtain the 
primary constraints defining the closed submanifold V = Im(J^£) ^ T*(J^7r), which are 

<^f^pU=^ ; ^f^{p\f-^,=^; (70) 

Let FCo'. J^TT — )• V. Then, the submanifold Wi = graphJ^£o = graphJ^£ is defined by 
Wi = [weWo: aw) = ^oiw) = ^tiw) = (l^fiw) = cl)^^\w) = O} 

where = — J^C p\, i = p — p. 

Next, we compute the tangency condition for the vector field X^ G X(Wo), given locally by 
(j69p on the submanifold Wi ^ Wo ^ W, by checking if the following identities hold 

L(X)eo^|^^=0 ; L(X)ef|^^=0 (71) 



L(X)<aS°^ , = ; \.{X)4>V , = . (72) 



a(0) 
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As we have seen in Section [2221 equations (j71|) give us the Lagrangian equations for the vector 
field Xq. However, equations (f72|) do not hold, since 



and hence we obtain two first-generation secondary constraints 

<^S'^=rfgi = o , 4^)=p0pi = (73) 

that define a new submanifold W2 ^ Wi. Now, by checking the tangency of the vector field Xq 
to this new submanifold, we obtain 

LTO0;^)=L(Xo)(pOgi) = O , L(Xo)4^)=LTO(p0p.') = -(p0)2, 

and a second-generation secondary constraint appears, 

4^ = {V\f = , (74) 

which defines a new submanifold W3 W2. Finally, the tangency of the vector field Xq on this 
submanifold gives no new constraints, since 

L(X„)</.(2)^L(Xo)((pO)2) = 0. 

So we have two primary constraints (|7Up . two first-generation secondary constraints (|73|) . and 
a single second-generation secondary constraint (|74p . Notice that these five constraints only 
depend on qf, and p\, and so they are ^2"Projectable. 

Notice that we still have to check (|7ip . As we have seen in Section 13.2.21 we obtain the 
following equations 

(-'^ - (-^)) (^^^ - {£) - ° 

(75) 



dq2dq. 



2 



Since we have already required the vector field X^ to be a semispray of type 1 in Wo, equations 
r6|) are satisfied identically and equations ([75]) become 



A long calculation shows that this equation is compatible if, and only if, 7^-^ = 0, for 1 ^ ^ ^ n 



dV 

That is, we have n first-generation secondary constraints arising from the tangency condition of 
Xq along Wi, thus defining a new submanifold W4 ^ W3 with constraint functions 

Observe that, since V is a function that depends only on t and g^, these new constraints also 
depend only on the coordinates t and q^., and thus they are /Sg-projectable. From a physical 
viewpoint, these constraints mean that the dynamics of the particle can take place on every level 
set of the potential with respect to the position coordinates. 
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Finally, we recover the Lagrangian and Hamiltonian dynamics from the unified formalism. 
For the Lagrangian solutions, using Proposition[71 we know that from the holonomic section ipo S 
r(p^) we can recover a holonomic section ipc = Pi^i^o & ^(Tf^) solution to equation (i36]l . In par- 
ticular, if V'o(t) = {t,qi{t),q{{t),qi{t),qi{t),p1{t),p\{t)), then i^c{t) = {t,ql{t),q\{t),ql{t),q\{t)) 
is a holonomic section solution to equations (j6ip . Now, bearing in mind the local expression 
(j63]) of the extendend Legendre-Ostrogradsky map, equations ([6T]) give the last n equations of 
the holonomy condition for '(/'£, which are identically satisfied since the holonomy condition has 
been already required, and the classical higher-order Euler-Lagrange equations 



dL 



dq^ 



A 



d dL 



dt dq{ 



A 



d2 dL 

+ 



dt2 dq^ 



= 0. 

■4>c 



For the Lagrangian vector field, from Lemma [3] and Theorem [51 we can recover from the 
semispray of type 1 Xo G X(>Vo) a semispray of type 1, Xc G X( J'^vr), solution to equations (|39|) 
(with Mf = plCWi)), and it is locally given by 

where are the solutions of equations (j77|) . 

One can check that, if the semispray condition is not required at the beginning and we 
perform all this procedure with the vector field given by (j68p . the final result is the same. This 
means that, in this case, the semispray condition does not give any additional constraint. 

Now, since C is an almost-regular Lagrangian density, for the Hamiltonian dynamics we must 
use the results stated in Section 15.31 and recover the Hamiltonian solutions passing through 
the Lagrangian formalism. For the Hamiltonian sections, by Proposition [T2l from a section 
ipo S r(/9|) solution to equation ([TT]) . we can recover a section ijjh = o p° o ip^ £ ^(t'o) 
solution to the equation (|l5]) . 



For the Hamiltonian vector fields, we know that there are semisprays of type 1 Xc S j£(J^7r), 
solutions to equations (|39|) . which are J^i2o-projectable on P4 = /52(W4), tangent to P4 and 
solutions to the Hamilton equation. 



7 Conclusions and outlook 



The objective of this work is to develop a complete and detailed geometric framework for describ- 
ing the Lagrangian and Hamiltonian formalisms of higher-order non-autonomous mechanical 
systems, and to give some applications of it. 

Our approach to the problem consists in extending the Lagrangian- Hamiltonian unified for- 
malism of Skinner and Rusk to this case, starting from the generalization of this formalism previ- 
ously made for first-order non-autonomous dynamical systems [5J and higher-order autonomous 
mechanical systems [39j. This enables us to derive the Lagrangian and the Hamiltonian for- 
malisms for these kinds of systems (in a natural way). We pay special attention to describing 
the equations of motion in several equivalent ways, using sections and vector fields in the bundles 
that constitute the phase spaces of these systems, and showing how the equivalence between the 
Lagrangian and the Hamiltonian formalisms is stated through the Legendre-Ostrogradsky map, 
which is also obtained in a natural way from the unified formalism. Our analysis is performed 
both for regular and singular systems. 
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As applications of our formalism, we study two physical examples: a regular system de- 
scribing the shape of a deformed elastic cylindrical beam with fixed ends, and a singular system 
describing a relativistic particle subjected to a generic potential depending on time and positions. 

The background geometrical tools that we use are higher-order jet bundles, in general, rather 
than the simpler and more usual trivial bundles (this particular case is also analyzed in the work) , 
since our aim is for this geometric framework to serve as a guideline towards a geometric model 
for higher-order field theories, which is free of the ambiguities present in their standard geo- 
metrical descriptions (concerning the definition of the Poincare-Cartan forms and the Legendre 
transformation). Some advances on this subject have been already obtained \10\ I46j . and we 
trust that our future work will contribute to completing them. 



A A particular situation: trivial bundles 

Assume that the bundle E M is trivial; that is, E = x Q, where Q is a n-dimensional 
manifold. In this case, we have that j'^vr = M x T^Q, where T'^Q is the kth order tangent bundle 
of Q (see [22] for details). The natural coordinates in this case are defined in the same way as 
in the general case, and are denoted by (t, Qq, qf, . . . , q^). In this case, the bundles involved in 
the construction are 



X T 



2k-l 



Q 



T*(J' 



X T*(T^-iQ) 



J TT = 



X T*(T'=-^Q) 



Thus, the higher-order restricted jet-momentum bundle is 

r = J TT Xjfe-i^ J TT 



X T— ^Q) Xkxt^-iq (M X T*(T^-iQ)) 



X T^'^-^Q XTfe-inT*(T'=-^Q) 



k-i, 



X Wa 



where Wa = T^^^-^Q T*(T'^ ^Q) denotes the unified phase space in the autonomous 

formalism. Natural coordinates in this bundle are the same as in the non-autonomous case, that 

,l2k-i'P%PA, ■■■ ■:Pa 



Remark: As we will see, in this particular situation the extended jet-momentum bundle is 
not needed. Thus, we denote Wr simply by W in this section. The differential forms Qr and fir 
(or, equivalently, Qq and ^lo) are also denoted by G and Q, respectively. 

We have the following diagram 



Pw, 




T*(T'=-iQ) 



where all the maps are the natural projections (see [39j for details). In coordinates, we have 



pM.{t, qf, qf,p\) = t , pw„(i, qf, qf,PA) = (qf, qf,p\) 
wi {qf, qf,PA) = (qf, qf) , p^2 (qf, qf,PA) = (.qf,PA) 
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Now we see how to construct the canonical structures in W, described previously, from the 
canonical structures in Wq. Let 9a E l7^(>Va) and Qa G ^^(VVa) be the canonical forms on Wa 
defined as 

where 9k~i and uJk~i are the canonical 1 and 2 forms on the cotangent bundle T*{T^~^Q). 

As stated before, the dynamics of the system is described by a Lagrangian density C G 
n^(R X T'^Q), with associated Lagrangian function L e C°°(M x T''Q). Then, if C G C°°(>Va) 
is the coupling function in the autonomous formalism ^39j . we can construct a globally defined 
Hamiltonian function in the following way: 

H = pl^C-L. 

Then, the forms Q E 17^ (W) and Q £ i7^(W) can be constructed as follows 
Q = PWa^a- Hp^r] , n = -de = pl^^Qa + dH A pI^t] . 

In local coordinates, bearing in mind the local expressions of co^-i and C: 

Ok-i = p\dqf , ujk-i = dqf Adp\ , C= p\qt+i , 

we have that the local expression for the forms and $7 are 

Q=p\dqt-{p\qti-L)dt , ^ = dqf A dp\ + d{p\qt+^ - L) A dt ■ 

that is, we obtain the local expressions given in (jlOp . 

The dynamical equations for sections and vector fields are now stated as in Section [31 and 
the local expressions are the same. There is only one difference: in Proposition [3] a connection 
in W is needed in order to split the presymplectic form O into the sum of a 2-form with the 
wedge product of two 1-forms (which are the differential of the local Hamiltonian function, and 
the volume form in R). In this case, we do not need to use such a connection, since the bundles 
are trivial and this splitting is natural. 
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